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Abstract 

We study exclusive vector meson photoproduction, yp —> V + p with V = J/'ij; or T, 
at NLO in collinear factorisation, in order to examine what may be learnt about the gluon 
distribution at very low x. We examine the factorisation scale dependence of the predictions. 
We argue that, using knowledge of the NLO corrections, terms enhanced by a large ln(l/£) 
can be reabsorbed in the LO part by a choice of the factorisation scale. (In these exclusive 
processes £ takes the role of Bjorken-x.) Then, the scale dependence coming from the remaining 
NLO contributions has no ln(l/£) enhancements. As a result, we find that predictions for the 
amplitude of T production are stable to within about ±15%. This will allow data for the 
exclusive process pp —* pTp at the LHC, particularly from LHCb, to be included in global 
parton analyses to constrain the gluon PDF down to x ~ 10~ 5 . Moreover, the study of 
exclusive J/U photoproduction indicates that the gluon density found in the recent global 
PDF analyses is too small at low x and low scales. 
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1 Introduction 


The present global PDF analyses (e.g. NNPDF3.0 [T), MMHT2014 [2], CT14 [3]) find that there is 
a large uncertainty in the low x behaviour of the gluon distribution. There is a lack of very low x 
data, particularly at low scales. Moreover, the gluon is determined at low x mainly by the DGLAP 
evolution of deep inelastic scattering, that is, not from a direct measurement^ but rather from 
the derivative di^/dlnQ * 2 . As a result, at low scales, the uncertainty on the gluon PDF is large 
for x 10~ 3 . On the other hand, the HERA data on diffractive vector meson photoproduction 
[11 El El El IS EH GUI HH H21021HH H5], 7 P —» V + p and the LHC data on the exclusive processes 
pp — >■ pV p [16, HJl IS] where V = J/ij) or T, and pPb —» pV Pb [19], sample directly the gluon 
distribution down to x ~ 10 -5 . For a review of exclusive vector-meson production at the LHC see, 
for example, Ref. [20]. 

However, J/-0 and T data are not used in the global PDF analyses. One reason is that the 
corresponding cross sections are described by non-diagonal (skewed) analogues of the PDFs, namely 
generalised parton distributions (GPDs), due to the different masses of the incoming photon and 
the outgoing V meson m The GPDs can be related to the PDFs via the Shuvaev transform in 
the low £ region [|22] j^] In this work, following [23], we make the physically motivated assumption 
that the input distribution has no singularities in the right-half Mellin-A r plane which implies that 
this relation holds at NLO to accuracy (9(£). Another reason is the dependence of the theoretical 
predictions on the choice of the factorisation scale, p, F . This problem has two independent parts. 
One part has a technical nature and the other is more physical. Let us discuss them in turn. First, 
the ‘technical’ problem, which is related to the convergence of the perturbative expansion at low 
£ and low scales. A good illustration is that the NLO amplitude for the exclusive high-energy 
7 p —* J/ij} + p process [23 J was shown to yield a cross section which varies by up to an order of 
magnitude for a reasonable variation of p F . This problem was emphasised recently by Wagner et 
al. [25] . The strong scale dependence arises because in the DGLAP evolution of low £ GPDs the 
probability of emitting a new gluon is strongly enhanced by the large value of ln(l/£). Indeed, the 
mean number of gluons in the interval A In fi F is [26] 

(n) ~ ^ hi(l/£) A In p? F , (1) 

7r 

leading to a value of (n) up to about 8, for the case ln(l/£) ~ 8 with the usual p F scale variation 
interval from p F /2 to 2/jl f . In contrast, the NLO coefficient function allows for the emission of only 
one gluon. Therefore we cannot expect compensation between the contributions coming from the 
GPD and the coefficient function as we vary the scale fi F . (At large £ the compensation is much 

Hn principle, measurements of Fj, would provide a direct determination, but the data are poor and do not reach 
low x. Moreover, the convergence of the perturbative series for is relatively poor. 

2 For GPDs, £ plays the role of the PDF variable x and is given by £ = ( p + — p' + )/[p + + p ,+ ), where p + and p' + 

are the light-cone plus-momenta of the in- and outgoing protons. 
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more complete and provides reasonable stability of the predictions to variations of the scale /ip.) 
In Section [2j we use the NLO contribution to fix a choice of the factorisation scale for the LO part 
of the amplitude, which allows the ln(l/£) corrections to be resummed. 

The second or ‘physical’ reason why there is a strong scale dependence of J/if and T photopro¬ 
duction is due to ‘defects’ in the presently available PDF sets which are obtained in the DGLAP 
global analyses of deep inelastic and related hard scattering data. That is, our studies indicate that 
there is a strong scale dependence in vector-meson photoproduction caused by the unexpectedly 
small LO contribution in comparison with the NLO correction. This occurs due to the very low 
input gluon density in the low x domain obtained in the recent global PDF analyses. We conclude 
that at the input scale, the low x gluon density is underestimated in comparison with that for 
quarks. The crucial observation is that the number of input gluons (which is parametrised freely in 
the global analyses) is found to be much less than the number of such gluons emitted by quarks - 
or, to be more precise, gluons associated with the quark (each quark carries a gluon held created 
by its colour charge). With the present global partons the NLO sea quark contribution to V meson 
photoproduction grows with decreasing x faster than the input gluons. As a result, at sufficiently 
high energy, the quark component of the NLO correction to V photoproduction approximately can¬ 
cels the main LO contribution which is due only to the gluon density. Therefore a small variation 
of the scale in the NLO component leads to a strong variation of the predicted cross section. 


Recently the LHCb Collaboration have presented data for open charm and beauty production 
in the forward direction HIII2E]. These data sample approximately the same kinematic domain as 
exclusive vector meson production; actually the x values are slightly larger, and the scale for beauty 
production is also larger. It was shown in j29| |30j that the data can be reproduced by NLO QCD 
using the present global PDFs. Note, however, that again there is a large sensitivity to the choice of 
factorization scale. Indeed, to reproduce the data one needs to take a scale approximately a factor 
of two larger than the natural choice /ip = \JQ 2 + 171 q • Since the optimal or appropriate scale for 
open heavy-quark production is not known at present, these data do not yet reliably determine the 
low x gluon at low scales. In the present paper we focus on exclusive vector-meson production. 
Moreover, our goal is not to present a new global PDF analysis or an explicit determination of the 
low x gluon, but rather to study the possibility to reduce the scale dependence of the prediction, 
to study the qualitative structure of the NLO amplitude and to give hints of the consequences for 
the behaviour of the low x gluon at low scales. 


In Section [2] we consider the NLO contribution originating from the light quark (singlet) GPD. 
We show that this part of the NLO amplitude allows us to choose a factorisation scale which sums 
the ln(l/£) enhanced contributions inside the GPD. The remaining part of the NLO contribution 
has no large ln(l/£) factors, and, as shown in Section 3.2, the compensation between the scale 
dependence of the GPDs and of the coefficient functions for T production makes the NLO result 
sufficiently stable for the data to be included in global parton analyses (just as for the predictions of 
processes at larger x). The study of the J/if photoproduction in Section 3.1 suggests that the gluon 
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density obtained in the existing PDF analyses is too small at low x and low scale. In Section [4] we 
present our conclusions. 


2 NLO corrections and the choice of fip 


Our aim in this section is to show how knowledge of the NLO contribution for yp —y V+p allows us to 
sum the ln(l/£) enhanced contributions by a particular choice of /ip in the LO part of the amplitude 
(and correspondingly in the NLO coefficient function), which, in turn, reduces the factorisation scale 
dependence of the predicted cross section. We outline the procedure in Section 2.2 , but first we 
comment on the NLO formalism. 


2.1 NLO formalism for ^p —> V + p 


The NLO formalism for yp — > Vp has been presented in [23]. However, before we can apply it to 
describe the high energy V photoproduction data, we must note one important correction. 

The LO partonic amplitude, Ag(x/£), for the subprocess y + (gg) —* V, can be calculated by 
considering the fusion of a photon with a pair of on-shell gluons with zero transverse momentum 
and physical, transverse, polarisations. If we use dimensional regularisation with D — 4 + 2e (as 
in [23]) to regularize the ultraviolet (UV) and infrared (IR) divergences, then the result is given, 
before dividing by the number of physical transverse polarisations of the incoming gluon, by 

A<°>(x/()= a,2(l +€). (2) 


In D dimensions, there are D — 2 = 2 + 2e possible transverse directions. Since the entire calculation 
must be performed in D dimensions, to properly average over the polarisation of the incoming gluon, 
we must divide ([2]) by this factor, and not simply by 2. Therefore for the LO partonic amplitude 
we obtain, contrary to [24] (prior to the first erratum), 

A£ ) (x/£) = a 8 . (3) 


The difference of a factor of (1 + e) does not alter the finite LO result (as e is taken to zero) but does 
lead to a different NLO result due to the differing counter-terms. Indeed, the extra factor (1 + e) 
presented in [23J generates extra terms ~ e/e in the counter-terms, A, giving them a form 


A = 


-+i+in(4 

e \p 2 


where 1/e = 1/e + y e — ln(47r) and y^ is the Euler-Mascheroni constant. However, with the LO 
partonic amplitude (|3]) we instead obtain counter-terms containing 
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It is precisely this alteration to the counter-terms which leads us to a different NLO finite result]^] 
This result agrees with that of EH with both errata. 

The details of our recalculation of the NLO result will be published separately [33J. Here we 
simply state the final result for the high-energy limit of the NLO part of the amplitude. In the 
high-energy approximation, W 2 My (where W is the c.m. energy of the incoming yp system and 
My the mass of the vector meson), the imaginary part of the amplitude dominates and the leading 
contribution to the NLO correction comes from the region £ <C |x| <C 1. The matrix element is 
given by 





(4) 


Here m is the mass of the c (for J/J>) or b (for T) quark and 

(5) 

( 6 ) 


Fg(x,€,fl) = V 1 - C 2 Hg(x,^n), 

Fs(x,^,fi) = y/l -£ 2 Hs(o;,£,/i), 


(for this unpolarised, forward process) with H g , Hs the gluon and quark singlet GPDs, respectively. 
The GPDs can be calculated from the diagonal PDFs, q(x,fi), q(x,y) and g(x,n), via the Shuvaev 
transforms: 
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V 0) 


4s(l — s) 
x + £(l-2s)’ 


3 This correction was also mentioned in Enng and discussed in more detail in |.'12j . 


(7) 

( 8 ) 

(9) 


( 10 ) 
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with an accuracy C?(£), see Ref. [2B]. Expressions ([8]) and (J9]) were used not only for the asymptotic 
NLO amplitude, Q, but also for the full expression for the complete NLO amplitude used in the 
numerics presented below. The LO coefficient function, 



OC s (^r) 


47ta/47 rae q (eye 7 ) f (Oi) 


NJi 


)v 


m° 


1/2 


( 11 ) 


where e* are polarisation vectors and (Oi)v is the NRQCD (non-relativistic QCD) matrix element 
of the cc —> J/-0 (or bb —>■ T) transition. In order to have a small relativistic correction to (Oi)v, 
we have to calculate the Feynman diagrams assuming that the charm/bottom-quark line has mass 
m = My/2 where V = T j34]. 


One can see directly from the high energy, leading ln(l/£), limit, given in Q, that for the choice 
Hf = m the asymptotic limit of the quark NLO contribution vanishes. It is this observation that 
will allow us (in Section 2.3) to claim that a suitable value for the factorisation scale in leading 
logarithm terms is /ip = m. 


Before we show how the NLO contribution allows us to fix the scale fip in the LO term, it is 
informative to first recall its structure and introduce the notation in terms of the LO prediction. 


2.2 General procedure 

We start with the LO contribution to yp —)• V + p. It is sketched in Fig. [lja) , and the amplitude 
is given by the convolution 

/ ! fj r 

—C< a >(x/()F'(x,(,p F ) = C<°> ® FM , (12) 

1 X 

where the sum over a = q, g is understood and Cq 0 ' 1 = 0 for this process. The coefficient function 
Cg 0 ' 1 is calculated using the non-relativistic vector meson wave function. In general, the relativistic 
corrections are not small for the J/-0 case. These corrections should be considered together with the 
three parton ( cc+g) component of the wave function; that is, accounting for the role of gluons which 
provide the interaction between charm quarks. For this process, it was shown by Hoodbhoy [21j that 
the consistent treatment of relativistic corrections may, to good accuracy, be effectively accounted 
for by choosing, in the non-relativistic formula, the charm quark mass m c = Mj/^/2^ After this, 
the remaining part of the correction is quite small (a few percent only). So, we may use the 
non-relativistic J/rp wave function with m c = Mj/^/ 2, to obtain a result with good accuracy. 

4 Strictly speaking, Hoodbhoy considered the electroproduction limit Q 2 M 2 This limit simplifies the 
calculation, but inspection of the proof indicates that the statement about a small relativistic correction should be 
reliable below this limit as well. 
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We return to discuss V = J/i/j or T. At LO the coefficient function Cg 0 ' 1 does not depend 
on fi F , whereas the low x gluon distribution F g depends strongly on the scale. In other words, 
based on ( 12 ), the exclusive V data measure g(x,fip), but we do not know the value of the scale 
at which it has been determined. How does this scale freedom arise? To obtain the LO result 
the coefficient function — the upper box in Fig. [jja) - is calculated with on-mass-shell gluons with 
transverse momenta l t = 0. In this collinear approach the factorisation scale /ip acts effectively as 
the ultraviolet (UV) cutoff of the logarithmic integral f oc In g 2 F over the gluon transverse 

momentum in the gluon loop of Fig. [IJa). Formally, in the LO collinear factorisation approach the 
value of /if is not known. In principle, the full result does not depend on fip since the higher-order, 
NLO, NNLO, ..., corrections compensate the effect of variations of fip- However, in reality the 
perturbative series is truncated and the compensation may not be sufficient to provide the scale 
stability of the theoretical prediction. 


On the other hand, we can go beyond the collinear logarithmic approximation by computing 
the l t integral accounting for the l t dependency of the hard 7 + gg -$-V matrix element, Ai, shown 
in the upper box of Fig. [jja). In comparison with the coefficient function (calculated with l t = 0) 
now the l t and l 2 dependence of Ai is included explicitly, and provides the UV convergence of the 
integral over l t . Formally, in collinear factorisation the difference between the pure logarithmic 
(ln/x|i) evaluation and the precise calculation of the gluon l t integral is treated as part of the 
NLO correction. This part of the NLO correction is of kinematic origin and is usually quite large. 
Fortunately, it can be moved into the LO component of the amplitude, noticeably reducing the 
remaining NLO correction. 


Instead of performing an independent calculation which accounts for the l t dependence, we can 
remain within the collinear approach and determine the value of the It integral given that the NLO 
coefficient function C^ of Fig. l](b) is known. Indeed, Fig. l](b) is the only diagram for the quark 
NLO coefficient function. In this approach the incoming quarks are assumed to be on-mass-shell 
and with zero transverse momenta but the loop integral over /, which contains the l t dependence, 
is calculated exactly. Since this is the same integral as that which occurs in Fig. [lja) we can use 
the result for C ^ to obtain a precise value, J, of the corresponding integral in the LO amplitude 
of Fig. [jja). After this we choose a scale hf — l^o which mimics the precise l t integration. That 
is, with a scale choice satisfying In/J q — J we have moved a large contribution from NLO to LO, 
and can continue to work in the conventional collinear approach, but now with a smaller NLO 
correction. 


5 A precise integration over l t was, in particular, implemented in 
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Figure 1: (a) The LO contribution to 7 p ^ V + p, showing the convolution (12). (b) The NLO quark 
contribution. For these graphs all permutations of the parton lines and coupling of the gluon lines to the 
heavy-quark pair are to be understood. Here P = (p + p ')/2 and l is the loop momentum. 


2.3 Transfer of part of the NLO to the LO contribution 


Since the above observation is crucial, let us demonstrate the procedure in more detail. At NLO 
level the LO+NLO amplitude at some factorisation scale pj may be expressed in the fornj^] 


.4<"7/) + ^l (1 7/) = 0 F( H ) + H F(^), 


(13) 


where the F s are the GPDs and where the coefficient function (7® does not depend on the factori¬ 
sation scale. Note that we are free to evaluate the LO contribution at a different scale p F , since 
the resulting effect can be compensated by changes in the NLO coefficient function, which then also 


becomes dependent on p F . Then eq. (13) becomes 


A^M + A^ip,) = cW®F(it F )+a,c£lM®F(n f ). 


(14) 


Note that although the first and second terms on the right hand side depend on pp , their sum does 
not (to 0(al)) and is equal to the full LO+NLO amplitude calculated at the factorisation scale pj. 


I 11 (13) the NLO coefficient function is calculated from Feynman diagrams which are 


in¬ 


dependent of the factorisation scale. How does the p F dependence of Crem in (14) actually arise? 
It occurs because we must subtract from the a s term which was already included in the LO 
contribution^] Since the LO contribution was calculated up to some scale p F the value of after 


6 For ease of understanding we omit the parton labels a = g, q on the quantities in (13) and the following equations. 
The matrix form of the equations is implied. 


''Simultaneously this subtraction also provides the infrared convergence of C^h 



















































the subtraction depends on the value /tp chosen for the LO component. The change of scale of the 
LO contribution from ///■ to /ip also means we have had to change the factorisation scale which 
enters the coefficient function from /if to /ip. The effect of this scale change is driven by the 
LO DGLAP evolution, which is given by 


A® (/ip) = 


C m - 111 | \ C <0 > ® V ) ® F(n f ), 


(15) 


where V denotes the skewed splitting functions. That is, by choosing to evaluate at scale /ip we 
have moved the part of the NLO (i.e. a s ) corrections given by the last term of (15) from the NLO 


to the LO part of the amplitude. In this way becomes the remaining /ip-dependent coefficient 
function CVem(/tp) of (14). In spite of the unusual form of (14), with two different scales /i/ and /ip, 
it is an exact equality at NLO and could, in principle, be generalised to higher orders]^] 


2.4 Large ln(l/£) terms resummed in the LO contribution 


The idea is to use the above procedure to reduce the scale dependence of the exclusive V pho¬ 
toproduction amplitude. Unfortunately the value of /ip is just one number, while C^ l \x /£) is a 
function of the ratio x/£. So we have no chance to nullify completely. Nevertheless we can 
nullify the most important NLO correction which is enhanced by the large value of ln(l/£). This 
contribution is generated by the integral J(dx/x) in the Crem(x/£, fXp)F(x, £, /i/) convolution of 
eq. (14) over the kinematic region of 1 x £. That is, we choose a scale /ip = /io which nullifies 
C^\x/^) iu the limit of x ^ £. It can be seen from Q that the scale /ip = /io = m ensures that the 
ln(l/£)-enhanced NLO corrections completely vanish for both the quark and the gluon components. 


From the NLO viewpoint the particular choice of /ip = m in the LO part is irrelevant. The 
corresponding order a s effect is exactly compensated by the remaining Orem term. However, the 
variation of /ip in the LO part affects not only the 0(a s ) terms but the higher-order a s contributions 
as well. Therefore, in this way, we resum the important large ln(l/£)-enhanced part of the higher- 
order a s corrections inside the parton distribution convoluted with the LO coefficient function and 
improve the convergence of the perturbative series]^] 


Actually our approach is rather close in spirit to the kt -factorisation method. Indeed, there, the 
value of the factorisation scale is driven by the structure of the k t = It (or the virtuality, Q 2 ) integral 


8 For example, if the NNLO contribution were known, then we will have three scales: /t/, /tp = /t LO , and /t NLO , 
where the NNLO correction to (14) takes the form a^Crem(/i LO ,/i NLO ) <8> F(nf). The scale /r NLO is fixed to nullify 


this term in the limit x £, and hence further reduce the sensitivity to variations of the scale [if. 

9 In particular, the most important factorisation scale dependence, enhanced by large ln(l/£), is caused by the 
double log terms, [a 8 In (/tp) ln(l/a;)] n , generated in the axial gauge by ladder-type diagrams. For GPDs this ladder 
contribution was studied in |36| . where it was shown, in the large ln(l/a;) limit, that the skewed splitting functions 
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in the diagrams of Fig. 00 In the /^-factorisation approach this k t integral is written explicitly, 
while the parton distribution unintegrated over k t is generated by the last step of the DGLAP 
evolution, similar to the prescription proposed in Refs. |37| .38]. Now, using the known NLO result, 
we account for the exact k t integration in the last cell adjacent to the LO hard matrix element. 
This hard matrix element A4, shown in the upper box of Fig. [l](a) , provides the convergence of the 
integral at large k t . In this way it puts an effective upper limit of the k t integral, which plays the 
role of an appropriate factorisation scale. 


The details of the prescription, for the case of the high-energy Drell-Yan process, were discussed 
in [39] • Indeed, Drell-Yan production of low-mass lepton pairs at high rapidity is another process 
for which the NLO prediction depends sensitively on the choice of the factorisation scale, unless 
the ln(l/a:) enhancements are first resummed in the incoming parton distributions. It was found in 
Ref. [39] that, after the scale fip = /i 0 is fixed for the LO contribution, the variation of the scale in 
the remaining NLO part does not noticeably change the predicted Drell-Yan cross section. In [39] 
it was shown how to calculate /io for the Drell-Yan process and, moreover, that the NLO prediction 
with /ip = /to is very close to the NNLO result. 


Returning to the earlier discussion in this subsection, it is clear from Q and (14) that the scale 
choice hf = tti = My/2 provides the appropriate resummation of the ln(l/£) enhanced terms, 
and, as a consequence, suppresses the remaining high-energy NLO corrections. In other words, the 
GPDs with = m (chosen in the LO part) include all the ln(l/£) enhanced contributions, while 
the remaining NLO corrections arise only from hard matrix elements corresponding to intermediate 
states of relatively small mass of about My. 


We can also explain this pictorially. As noted above, the fip dependence of the NLO correction 
is caused by the subtraction of the contribution generated by the evolution equation. In particular, 
the correction which (i) depends on /ip and (ii) is enhanced by ln(l/£), is generated by the ladder- 
type diagram^ shown in Fig. [2] The choice of the factorisation scale (ipi determines which part 
of the diagram is attributed to the evolution of the PDF and which to the hard matrix element. 
By choosing the value of fip = fio — My/2 we include all the ladder cells in the LO part, so that 

are proportional to 1/x. So these double log terms come from DGLAP integrals of the form 


(/•* f 



dkl_ 1 r k i d*?\ / r 1 dxn f 1 

kn- 1 J k'f J \J X X n J Xn 


dx n -1 

%n— 1 


r 1 dx\ 

lx 2 Xl 


[a s ln(g F ) ln(l/x)]"/(n!) 2 


where we have strong k 2 = k 2 and x ordering. Thanks to the strong k t and x ordering in these double log LO 
integrals, the correct upper limit, g F , in the first integral automatically provides the exact resummation of all the 
terms in the double log series. 

10 We stress again that, in the high energy (x £) contribution, the form of the integral over the gluon loop 
momentum l t is exactly the same in both the quark (Fig. 0b)) and gluon (Fig. 0 b) but with the quark lines 
replaced by gluons) NLO contributions. Therefore the scale g F = Uo simultaneously nullifies the high energy quark 
and gluon contributions. Note that this is only true after the corrections discussed in Section 12. II are included. 

11 Ladder diagrams occur in the axial gauge which is conventionally used to calculate the GPDs. 
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Figure 2: A Feynman diagram which illustrates different possibilities for the division of cells between the 
parton evolution and the hard matrix element, dependent on the choice of the factorisation scale, fipi- 


the remaining NLO matrix element will not receive any contribution from ladder-type diagrams (at 
small £). Note that the ladder diagrams would contribute to the NLO matrix element (from /ip to 
Po) if we were to choose a scale /ip < ho- For the choice /ip > /i 0 the NLO matrix element would 
contain the ladder diagrams with the negative sign in order to compensate the extra contributions 
(from /i 0 to /ip) included in the LO part. 

Since we consider an exclusive process, the mass of the final state is fixed, but the intermediate 
states, corresponding to the discontinuity which gives the imaginary part of the amplitude (depicted 
using a dotted line in Fig. [2]), depend on the choice of scale. For larger values of /ip the mass of 
the allowed intermediate states in the matrix element becomes smaller. Of course, it is impossible 
to absorb inside the evolution all intermediate states other than the heavy quark pair QQ state 
which produces the V : for example we cannot absorb the effects of the uppermost gluon in Fig. [2] 
emitted by the upper heavy quark. Nevertheless, when we choose the appropriate scale /ip = My j 2, 
the mass of the remaining intermediate states becomes of the order of My. This choice will move 
ln(l/£) contributions from the matrix element into the GPD. It will therefore provide a much more 
stable final result, since the remaining NLO contributions then cannot be enhanced by the large 
values of ln(l/£). 

It should be emphasised that the asymptotics of the NLO amplitude is used only to determine 
the effective scale /ip. In all our further numerics we use the full expression for the complete NLO 
amplitude as given in [21, 32j [33]. 
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3 Can exclusive vector meson production data be included 
in a global PDF fit? 


In this section we present our results for the exclusive 7 p ^ J/i/> + p and yp —> T + p processes as 
functions of the factorisation scale pj, and the renormalisation scale pp. The cross section relies 
sensitively on the gluon PDF at small x, in a kinematic regime where it is very poorly known, 
as well as depending of the choice of scales. For these quasi-elastic processes we present just the 
dominant imaginary part of the amplitude. The real part of the coefficient functions has been 
calculated exactly, both via a dispersion relation [23] and by directly computing the real part of the 
loop integrals [321l33|. it is non-zero in the time-like region |x| < £. Therefore, to compute the real 
part of the amplitude directly, we need also the GPDs in this region. Unfortunately, the Shuvaev 
transform is not valid in the time-like region [23]. Nevertheless, the real part of the amplitude can 
be included via dispersive methods on the level of the amplitude. However, giving a prediction 
of the full cross section is not our objective here. Rather our aim is to study the stability of the 
perturbative predictions and to investigate whether or not we can determine optimum scales so 
that experimental data for these processes, and the related pp pVp processes, can be included 
in global PDF (collinear) analyses to constrain the gluon PDF at low x, in a domain for which, at 
present, there are no data. For this goal it is sufficient to work with the more simple imaginary 
part of the amplitude. 


We recall the main result of Section 2. We start from the key equation, (14), that is 


+ AW( N ) = C<°> ® FM + a,CgLM ® F( N ) . (16) 


To obtain this result we had introduced by hand a new scale jip in the LO term and showed that 
the choice pp — m — My/ 2 allowed us to resum and transfer all the enhanced large ln(l/£) terms 
from the NLO contribution to the LO term, leaving a much smaller remaining NLO contribution 
ol s C'rem ( Pf ) <S> F{pf). There is still a pf scale dependence, but now this should be much weaker. 


3.1 The process 7 p J/t/j +p 

To obtain predictions for exclusive J/i/> photoproduction we are working at scales close to the input 
GPDs used in the calculation. The results for the scale variation of the imaginary part of the 
amplitude are shown in the two plots of Fig. [3] In each plot we show separately the LO and NLO 
contributions to the amplitude J* 2 ] The left plot shows how these contributions change if we vary all 
the scales pf = pp = pp = p simultaneously, taking p 2 = m 2 / 2, m 2 , 2m 2 . In the right plot we 

12 We choose to use CTEQ66 partons [40] since the gluon obtained is positive definite and since they were used 
in our earlier works on the subject BUGS]- Moreover, since we are emphasising the general procedure rather than 
making quantitative predictions, the choice of any particular, reasonable PDF set is not important. 
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W [GeV] W [GeV] 


Figure 3: Predictions of Im A/W 2 for y p —> J/if) + p as a function of the 7 p centre-of-mass energy W, 
produced using CTEQ6.6 partons [4U] with scales fj, = hf — /J>f — fJ>R (left panel) and p = Hf — hr 
with fi 2 F = 2.4 GeV 2 (right panel). The bottom (top) set of curves corresponds to the Born (1-loop) 
contribution with the scale variation p 2 = 1.2, 2.4,4.8 GeV 2 . The dot-dashed, solid and dashed lines 
correspond to the low, central and high values of the scale /j, respectively. Note that for the left panel, 
the bands overlap for energies bigger than about 70 GeV. 


fix /ip at the optimum scale fip — m and vary /if — fi R = fi. The result is dramatic. The transfer 
of the ln(l/£) terms from the NLO to the LO contribution has significantly reduced the py scale 
dependence. 

But there is another, more severe problem. The LO contribution is dominated by the NLO 
contribution of opposite sign. Thus the imaginary part of the quasi-elastic yp — >• J/'ijj+p amplitude 
changes sign when the NLO contribution is added. As it stands, this result is in contradiction with 
modelling the interaction as an elastic forward scattering where the imaginary part of the amplitude 
is positive. What is happening? The explanation is interesting. In general, the global DGLAP 
PDF analyses start from input forms which are completely arbitrary and, moreover, neglecting any 
constraints, know nothing about the structure of the evolution at low Q 2 . It is then found that the 
gluon PDF tends to be small, or even negative, in the low Q 2 , low x (10 -4 < x < 1CV 2 ) domain, see 
Fig.0 Clearly due to the lack of data constraints in this domain the gluon is not reliably known. 

Let us study the over-simplified case with an input gluon g(x) = 0 so that at the input scale 
we have only quark PDFs. For yp — y J/'ip + p the quark contribution only appears at NLO. The 
imaginary part of this NLO contribution is negative with respect to the normal (g 7 ^ 0) LO term. 
Indeed, when calculating the NLO coefficient function we must subtract the contribution which is 
already generated by LO evolution. However, the subtraction is performed purely formally in order 
to avoid the infrared divergence and does not account for the actual value of the gluon density used 
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Figure 4: The gluon distribution at Q 2 = 1.21 GeV 2 as determined by two recent global analyses [HE]. 
Figure produced using [[421 [45]. 



to calculate the LO term. Therefore, at low scales, where the LO term is suppressed by much too 
small ‘unphysical’ gluon PDFs obtained from the global fits, the NLO term generated by the quark 
is not just a correction, but is the dominant contribution. As a consequence, the imaginary part 
of the quasi-elastic amplitude becomes negative (in comparison with the LO contribution) and the 
observed behaviour reveals a lack of stability of the perturbative series at this order. Of course, 
at high scales where evolution, and not the input, determine the PDFs, we obtain a sensible NLO 
prediction. 

This simple consideration demonstrates that if yp —>■ J/ip + p and pp — > p + J/ip + p data 
could be included in a global fit, then they would put a strong constraint on the low x input gluon 
distribution. It is not simply that g(x) must be positive, but that actually the input gluons cannot 
be smaller than the density of gluons emitted by the quarks before the beginning of the evolution, 
when the parton virtuality Q < Qq. If J/ip data would have been included in the collincar global 
parton analyses, this constraint on the low x input gluon would have been automatically satisfied. 

Indeed, it is seen from the second plot of Fig. [3] that after the ln(l/£) enhanced corrections 
are resummed by fixing /ip = m, the stability of both the LO and the NLO components of the 
amplitude under the scale variations are much better. From this viewpoint the J/ip data may be 
included in a global PDF analysis. The only problem is that the gluon density obtained in the 
existing PDF analyses is too small at input scales and low x (10 -5 < x < 10~ 2 ), so that we get the 
wrong sign of the imaginary part of the amplitude. This means that if J/ip data were included in 
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the global PDF analyses, we must get a larger gluon density at low scales]^] Larger gluons in this 
kinematic domain will increase the LO component of the J/-0 photoproduction amplitude and will 
provide the correct sign for the whole LO+NLO amplitude. 

3.2 The processes ^p —> T ± p and pp -A pTp 

The top two plots of Fig. [5] show the results for the scale sensitivity of yp —* T + p, analogous to 
those of Fig. [3] for yp — > J/ip+p. Again we see that fixing the scale pf = Mr/2 reduces the pf scale 
uncertainty. The lower two plots show that part of this stability arises because the change caused 
by variation of p R is to some extent compensated by an ‘opposite’ change due to the variation of pf. 
If we were to choose a non-optimal scale, say for example, p 2 F = 2m 2 , then the scale variation turns 
out to be about twice as large as for the optimal choice p 2 F = m 2 . Due to the larger T mass we are 
now working at scales with more perturbative stability. The LO contribution is partly cancelled by 
the NLO term, but is not dominated by it. If we take the upper right or lower left plot of Fig.[5]then 
the pf scale uncertainty of the amplitude is about ±15% and ±25% respectively. As a result the 
data for exclusive T production can be used in global PDF analyses to probe the gluon distribution 
down to x ~ 10 , in the case of LHCb kinematics. The calculation of exclusive pp —>■ pTp is 

described in [35]. It is based on the sum of the two diagrams shown in Fig. [6] For an T produced at 
large rapidity y , the dominant contribution is from the diagram with the larger yp centre-of-mass 
energy W + , which depends on the gluon density at x ~ My e~ y j\fs. The small contribution from 
the other diagram, with much lower energy VF_, may be estimated from the existing HERA data. 

If we use the central values of the presently available PDF sets obtained from global analyses, 
then we find cross section predictions which are about a factor of four below the HERA exclusive T 
data mmm- This indicates that future PDF global analyses with T data included will, just as 
for J/i[j data, require a larger gluon distribution at low values of x. While the exclusive T process 
samples the PDFs at scales Q 2 ~ M\/ 4, the larger gluon required at these Q 2 will affect the gluon 
densities at all values of Q 2 via the evolution. 

3.3 Note on the alternative Ay-factorisation approach 

The collinear DGLAP global analyses tend to result in gluon distributions with valence-like x 
distributions at low scales, see Fig. [4j On the other hand, the approach used in the JMRT paper [35] . 
to study exclusive vector meson production, is free from this problem. There, the ‘NLO’ prediction 
was not calculated as a correction from NLO Feynman diagrams in collinear factorisation, but 

13 Most probably the pure DGLAP global analyses should take into account absorptive corrections, which are not 
negligible at low x and low Q 2 . The existing very small gluon densities in this kinematic domain are a way of 
mimicking these absorptive effects. 
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W [GeV] W [GeV] 

Figure 5: Predictions of Im A/W 2 for 'yp — » T + p as a function of the 7 p centre-of-mass energy W, 
produced using CTEQ6.6 partons [40j. The scales are /1 = n F — A 4 / = Hr (top left), /x = /j,f — /x R 
with fi 2 F = 22.4 GeV 2 (top right), /i = /// with /i 2 F = /i 2 R = 22.4 GeV 2 (bottom left), n = /j, R 
with /i 2 = /j,p = 22.4 GeV 2 (bottom right). The top (bottom) set of curves corresponds to the LO 
(LO+NLO) contribution with the scale variation /jl 2 = 11.9,22.4,44.7 GeV 2 . The dot-dashed, solid and 
dashed lines correspond to the low, central and high values of the scale /j, respectively. The large GPD 
uncertainties are not shown. 


16 


























Figure 6: The two diagrams describing exclusive T production at the LHC. Diagram (a), the W + 
component, is the dominant contribution for an T produced at large rapidity y. Thus data for pp —» pTp 
allow a probe of very low x values, x ~ Mr e _y /yfs; recall that in the (dominant) imaginary part of the 
Born amplitude we have x = 


approximated by taking the full integral over the gluon kt in Fig. [lja) , including an ansatz for the 
now kt dependent gluon. The convergence of this explicit kt integral provides effectively the ‘optimal’ 
value of the factorisation scale and in this way sums all ln(l/£) enhanced terms, accounting for a 
large part of the NLO corrections. It does, of course, not include corrections from diagrams which 
do not have a ladder structure, however, the contribution from the light quarks, Fig. [ljb) , is already 
included in the unintegrated gluon distribution used in JMRT [35] . Clearly, then the quasi-elastic 
amplitude is positive definite. Let us briefly describe how our ‘NLO’ predictions for yp —)■ T + p 
data were made. First the incoming gluon distribution was fitted, within this approach and using an 
ansatz for the gluon based on the important double-logarithmic dependence on ln(l/a;) and In Q 2 , 
to reproduce the J/-0 data from HERA and the LHC. Then we proceeded to T production using 
our fitted gluon. We verified that, in the relevant kinematic domain, this gluon reproduces the NLO 
DGLAP evolution to good accuracy. Therefore it was not surprising that the ‘NLO’ predictions for 
the T cross section, as a function of W, agreed well with the LHCb data [18j when they became 
available. 


4 Conclusions 

Here we have been concerned about the inclusion of exclusive vector meson production data in 
global PDF analyses in order to probe the gluon density at small values of x; that is, in the x < 10~ 4 
domain. Note that in this domain the input gluon PDF is freely parametrized, and is found to have 
a tendency to be valence-like with large uncertainties. It was hoped that the situation would be 
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changed as data for exclusive vector meson production, V — J/ip and T, which are very sensitive 
to the gluon at small x, became available at HERA and the LHC; particularly measurements of 
the exclusive process pp —> p V p at the LHC with V detected at large rapidity. Why has this not 
happened (so far)? 

For J/-0, the main problem is, at first sight, the very poor convergence of the perturbative 
expansion in the collinear approach. This can be seen, e.g., from Fig. [3j where the NLO contribution 
is comparable to the LO term and opposite in sign. However, we argued that this large NLO 
contribution, in comparison with that of the LO, reflects not the poor convergence, but rather that 
the global PDF analyses find a gluon density which is too small at low x and low scales. For this 
reason we find a LO contribution to exclusive J/ip photoproduction amplitude which is too small. 
We argued that the input gluon in the collinear approach, used in the global PDF analyses, should 
not be parametrized freely, but should be subject to some constraints. We noted that working in 
the physical, Ay-factorisation type scheme would avoid these problems. 

For exclusive T production the situation is much better. We found that the optimum factori¬ 
sation scale is much higher, the perturbative expansion at NLO in collinear factorisation converges 
well, and the remaining mild scale dependence of the predictions (~ ±15% on amplitude level) 
means that data for pp —y pTp can now be included in the global PDF fits to determine the gluon 
in the low x regime for the Erst time. 


It is appropriate to list the theoretical uncertainties of the present calculation. Although the 
leading double log terms have been resummed correctly to all orders, there still exist remaining 
NNLO and higher contributions, which are unknown at present. When the full NNLO amplitude 
is known we showed how the scale uncertainty can be further reduced. Next, we consider the 
accuracy of the expressions which relate GPDs to conventional PDFs, (J8]) and ([9]). These relations 
are based on conformal invariance and the equality of the Gegenbauer moments of GPDs to the 
Mcllin moments of PDFs. Due to the polynominal property [21] the accuracy of this equality of the 


moments is (9(£ 2 ), providing an O(g) accuracy for ([8]) and (Jt)]), see [23]. However, recall that the 
Shuvaev transform assumes the absence of the singularities in the right-half Mellin-iV plane for the 
input distribution. This assumption is physically reasonable since in the Regge approach there are 
no singularities in the right-half (j > 1) plane in the space-like (|x| > £) domain where we actually 
work. Nevertheless, whenever possible, this assumption should be checked. One check is that the 
predictions for the GPD/PDF ratio obtained in this way are in a good agreement (especially for low 
x gluons, which are the most important for exclusive J/^ or T production) with the NLO results 
of |44j coming from a fit to deeply virtual Compton scattering (DVCS) HERA data. Finally, the 
relativistic correction to the vector meson wave function, which is discussed in Section 2.2[ is not 
expected to be significant pkTj . 

We noted that an alternative probe of the low x gluon at low scales has been considered in 
There they study the data for charm (and beauty) production obtained in the forward 
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direction by the LHCb Collaboration [23 US\- Again it is found that the predictions strongly 
depend on the choice of factorization scale. Their prediction is lower than the data if the natural 
scale fj/p = m 2 Q +p 2 is chosen, and it is found that a larger pp is needed to reproduce the data using 
the existing global PDFs. Thus our expectation of a larger gluon at low x is not in contradiction 
with the LHCb charm (and beauty) forward data. 
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